Shot-noise-driven escape in hysteretic Josephson junctions 
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We have measured the influence of shot noise on hysteretic Josephson junctions initially in macro- 
scopic quantum tunnelling (MQT) regime. Escape threshold current into the resistive state decreases 
monotonically with increasing average current through the scattering conductor, which is another 
tunnel junction. Escape is predominantly determined by excitation due to the wide-band shot noise. 
This process is equivalent to thermal activation (TA) over the barrier at temperatures up to about 
four times above the critical temperature of the superconductor. The presented TA model is in 
excellent agreement with the experimental results. 
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Shot noise and full counting statistics (FCS) in meso- 
scopic conductors are currently intensively studied, be- 
cause such measurements yield fingerprints of conduc- 
tion mechanisms of these scatterers y, La |j] . Yet di- 
rect measurements of the second and higher moments of 
current or voltage are typically much more difficult and 
prone to errors as compared to the measurement of DC 
transport properties. This is because, in order to extract 
information about fluctuations, one generally needs to 
perform high frequency measurements on remotely con- 
nected samples at cryogenic temperatures, which is in 
practise a non-trivial task. Therefore, detectors of noise 
and FCS that operate directly on-chip near the noise 
source are of great importance. A recently proposed 
read-out device of these fluctuations is a Josephson junc- 
tion (J J) threshold detector 0,0. A mesoscopic non- 
hysteretic Josephson junction in the Coulomb blockade 
regime is another choice 0, H, • 

The dynamics of the JJ can be described as that of a 
phase (ip) particle in a tilted cosine potential, see Fig. ^ 
Under the influence of equilibrium environment fluctu- 
ations at low temperatures the particle resides in the 
quantum mechanical ground state, from where the es- 
cape mechanism is tunnelling (MQT) through the bar- 
rier. At higher temperatures it assumes a nearly thermal 
population of the states, and the escape mechanism is 
predominantly thermal activation (TA) over the barrier 
top. Below we discuss the regime where excited states 
are not accessed by the influence of the thermal noise 
of the dissipative environment, which is an appropriate 
assumption based on our present measurements in the 
absence of shot noise and those of Ref. 01 m the same 
set-up at low temperatures. 

The influence of noisy current on escape characteris- 
tics is typically considered in two different limits: (i) 
adiabatic (quasistationary) regime, where fluctuations of 
bias current that change the tilt of the cosine potential 
are so slow and weak that the phase particle remains in 
the ground state of the metastable well, and the escape 
rate is dictated by MQT from this varying ground state 
BEL or (h) (nearly) resonance excitation limit, where 
most of the escape events are consequences of shot-noise- 
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FIG. 1: (a) Measurement configuration. A split JJ in MQT 
regime measures the current fluctuations due to the shot noise 
of the scattering junction in the right arm. (b) Processes lead- 
ing to escape from zero-voltage (supercurrent) state into resis- 
tive state, (c) Sketch of current pulses Ip, with superimposed 
noise of 7n, and (d) the associated escape characteristics. The 
central quantities relating to the pulses and escape histograms 
have been indicated in the figure. 



driven excitations of the phase particle into higher levels 
in the well. Due to the high frequency resonance char- 
acter of our on-chip circuit, it turns out that the latter 
approach is more appropriate in describing the dynamics 
of the detector in this case. 

In this Letter we show that a hysteretic Josephson 
junction in MQT regime can measure high frequency 
current fluctuations generated by a mesoscopic scatterer 
in non-equilibrium. We present a model, in which the 
noisy current excites the phase particle in a metastable, 
nearly parabolic well out from the ground state, and sub- 
sequently the particle escapes from the well. This occurs 
via thermal activation at an equivalent temperature of 
the phase particle which is determined by the competi- 
tion between shot-noise excitation, and relaxation due to 
the dissipative environment. This model is in excellent 
agreement with our experimental results, and we are able 
to investigate Josephson phase dynamics up to equivalent 
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temperatures which are about four times higher than the 
critical temperature of the superconductor. This is the 
temperature responsible for phase dynamics: the elec- 
trons and lattice of the superconductor are only weakly 
coupled to this subsystem, and the former two remain 
essentially at their base temperature. 

Our measurement configuration with an electron mi- 
crograph of a split JJ detector is shown in Fig. ^ 111 a 
typical experiment two currents, In and ip, are injected 
into the circuit through large resistances (> 1 Mil) at 
room temperature. Current In is applied constantly and 
it runs through another tunnel junction at a distance of 
120 /im from the detector. This junction plays the role of 
the shot-noise source in the circuit. The DC-component 
of this current, In, returns through a long (3 mm) and 
narrow (width 2 fim) superconducting line such that no 
DC current due to In passes through the detector. We 
verify this balance frequently during the measurements. 
Thus, only the fluctuations of 7n are admitted through 
the detector. We probe these fluctuations by applying re- 
peatedly trapezoidal current pulses of height Ip through 
the detector as shown in Fig. |Tj (c). These pulses are 
slow, typically 10 /is - 10 ms long and with long transient 
times at the beginning and at the end of the pulse to se- 
cure adiabatic response to them. Typically 1000 pulses 
at each value of Ip are repeated, and the escape prob- 
ability P{Ip) is obtained as the fraction of those pulses 
that lead to escape from the supercurrent state. 

Let us discuss the dynamics of a hysteretic JJ in more 
detail. In the common RCSJ-model (Resistively and Ca- 
pacitively Shunted Junction model), the phase particle in 
the tilted cosine potential is cither trapped and oscillat- 
ing in one of the wells, or alternatively it runs down the 
potential. The latter regime is called resistive state. The 
measurable quantities, i.e., current / and voltage V, are 
related to ip via the Josephson equations: / = Icsmip 
and 2eV = fo^r, where Iq — j^Ej is the critical cur- 
rent of the junction and Ej is its Josephson coupling 
energy. The quality factor Q = u p RC is assumed to be 
Q 3> 1 to assure hysteretic dynamics. This means that 
once the particle leaves the well, it runs freely such that 
the voltage across the junction is close to twice the en- 
ergy gap of the superconductor, and it can be retrapped 
to the zero-voltage supercurrent state only when cur- 
rent is lowered virtually to zero. Here R is the resis- 
tive shunt of the junction, C is the parallel capacitance, 
and w p = \/8EjEcqo/h is the plasma frequency, where 
Eq = e 2 /2C and qo = -\/2(l — I/Ic)- Quantum mechan- 
ically the phase particle has energy states in the nearly 
parabolic wells, with energies close to those of a harmonic 
oscillator: E n ~ (n + l/2)frui p . 

Due to the nearly harmonic potential, transitions be- 
tween the neighbouring levels are dominating. The rates 
of these transitions, see Fig. ^ (b), are determined by 
the spectral density 5/(w) = /^(/(^/(O)) 
of the current noise at the corresponding level separa- 
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and 7j-i.j — (i/2hu)j j-iC]Si(+u)j for excitation 
and relaxation, respectively [21 |l3j . The noise spectrum 
has two contributions, one due to the equilibrium envi- 
ronment, S™ v , and the other due to shot noise, S} hot , 
which add incoherently: 57 = S c j nv + S} hot . At low 
temperature, fc B T -C huj p , the thermal excitation is 
strongly suppressed, SJ nv (— ojjj-i) — > 0, and the relax- 
ation is due to zero point fluctuations, Sj nv (+uijj-i) — > 
2Hu)jj—iRjeY(u)jj—i), where Y(uS) is the admittance of 
the circuit surrounding the Josephson junction. The shot 
noise responsible for transitions in the Josephson junc- 
tion can be described as Sj hot (±ujjj-i) = Felys. Here 
In is the average current through the scatterer, and F 
is the "Fano factor" of the noise source and the circuit 
surrounding the J J, at the frequency corresponding to 
level separation. Factor F can be calculated for a known 
experimental circuit. Here we determine it experimen- 
tally as a fit parameter. It would be unity in the case of 
Poissonian tunnel junction source [lj and if all the noise 
current would run through the detector junction. Com- 
bining the results above, we have 
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The level dynamics described by Eqs. and @ 
can be described by the equivalent temperature 
T* and effective quality factor Q* by requesting 
7j,j_l = J ^Q' 1 [coth(7kjj,j_i/2fc B r*) - 1] and Tj-ij = 
j^i[coth(^ J _i/2fc B r*) + 1]. This yields Q* = Q 
and, with Ujj-i — w p 
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Following the standard results of the decay from a cu- 
bic mctastable well, we can infer that thermal activa- 
tion is the dominant escape mechanism provided T* > 
To = hujp/2i:k^. This yields a condition for the valid- 
ity of the TA model: arcoth(l + QFeln/h^lC)/^ < 1. 
This is fulfilled at all currents Jn that were employed in 
the experiments, as we will show below. With this pro- 
cedure it is then straightforward to obtain the switch- 
ing probability of the threshold detector, in the limit 
of many levels in the well, as -P(-Tp) = 1 — exp(— TAt). 
Here T = ^ exp(-AC//fc B T*) is the standard TA es- 
cape rate, At is the length of the current pulse, and 
AU = AV2Ej(l - 7p// c ) 3/2 /3 is the height of the po- 
tential barrier at the particular bias point Ip. The ex- 
pressions above allow us to evaluate the position Jo. 5 an d 
the width AI = Jo g — Iq.i of the escape threshold, where 
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TABLE I: Parameters of the two samples. 

Sample I c C I c C 

detector detector scatterer scatterer 

A 3.9 pA 230 fF 600 nA 40 fF 

B 1.5 M 230 fF 15 nA (1) 10 fF (1) 

90 nA (2) 40 fF (2) 



I x is defined by P{I X ) = x, see Fig. n](d). These approx- 
imate results were compared to those obtained by full 
level dynamics calculation [12J numerically. They were 
identical, although there the results depended sensitively 
on the exact values of the level positions and widths near 
the barrier top. 

In the limit of large noise currents, Jn 3> we 

obtain T* ~ S FeIf U , which can be written in a more 

familiar looking form T* ~ 2fcs AR ' ^his ^ s a vau d a P~ 
proximation in a wide range of experimental parameters. 
The expressions of TA escape and the current dependence 
of the barrier height combined with this approximate ex- 
pression of T* finally allow us to approximate Jo. 5 and 
AJ as 

h,/ic - i - (^mH^r 3 (QF^r /3 

(4) 

and, further assuming UpAt 3> 1 yields 

lnK 9 ) - ln( Ko . x ) Eg jj_ )a /3 (5) 
' C [i21n(^)]i/3 ^ Ej eu p J 1 ; 

Here k x = — ln(l — x). 

We report on data of two samples, A and B, which 
were fabricated by standard electron beam lithography 
and shadow evaporation with aluminium as the super- 
conductor. The samples were measured via adequately 
filtered signal lines in dilution refrigerators at bath tem- 
peratures of 30 mK - IK. Sample A had a split JJ detec- 
tor and one scattering junction as in the scheme of Fig. 
n Sample B had a single JJ as a detector, and it had two 
scattering junctions located symmetrically with respect 
to the detector (at a distance of 120 /jm), and with re- 
spect to two long injection lines. In this sample the two 
scattering junctions were made intentionally very differ- 
ent to check the invariance of the results with respect to 
junction properties. The parameters of the two samples 
are listed in Table HJ 

Data in Fig. [21 (a) show results of a control experiment 
on Sample B at the base temperature T ~ 30 mK, where 
the trapezoidal pulse current J was injected through (i) 
one of the long injection lines, (ii) through scatterer 1, 
and (hi) through scatterer 2. The histogram of case 
(i) lies at higher currents than those of (ii) and (iii), 
which in turn overlap practically with each other. This 
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FIG. 2: Histogram positions under shot noise injection, (a) 
Escape histograms of Sample B when 800 ^ s long current 
pulses have been injected through the long injection line (open 
circles), and through the two noise sources (1 - squares, 2 - 
filled circles), respectively. The lines are the corresponding 
theoretical results, (b) Experimental results (circles) on the 
switching threshold current 7o.s against the average current 
In through the scatterer 2 in Sample B. Pulse length was 
At = 800 fis. The solid line is again the result of the used 
theoretical model. The inset shows the IV curve of the noise 
source, (c) Similar data as in (b) but for Sample A. The 
rising curve shows the corresponding equivalent temperature 
T* , which ranges from 2 to 5 K in this measurement. 



demonstrates that shot noise tends to push the thresh- 
old towards lower values of current, as predicted, e.g., 
by Eq. (0J. Furthermore these data demonstrate that 
the noise is predominantly, and equally in (ii) and (iii), 
generated by the scattering tunnel junctions. The calcu- 
lated lines run through the corresponding experimental 
histograms assuming MQT without any noise in (i) and 
TA with T* evaluated from Eq. © with QF = 5 in 
(ii) and (iii) . The latter value is realistic in terms of the 
expected O ~ 10 and F < 1. 

Figure [21 (b) shows data on Sample B at the tem- 
perature T = 40 mK. The threshold current Jo. 5, i-e., 
the pulse current at which the switching probability is 
50%, has been plotted as a function of Jn through the 
scatterer (2). The data, shown by solid symbols, follow 
the prediction of the model presented above, again with 
QF = 5. The initial plateau in the experimental data 
below Jn — 100 nA arises because the scattering junc- 
tion is not in the linear quasiparticle tunnelling regime 
here. The inset of Fig. [21 (b) shows the IV curve of this 
junction. 

In Fig. [21(c) we plot, in addition to similar data as in 
(b), the equivalent temperature T* in the measurement 
of Sample A when current Jn is varied between 2 fiA 
and 8 [iA. There are two interesting points to note, (i) 
We are able to study Josephson dynamics of the JJ up 
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FIG. 3: The width of the histograms when shot noise and 
thermal noise are applied, (a) The width of the grey zone, 
AI, against the average current Jn through the scatterer. The 
dots are the experimental results and the solid line is the 
result of the theoretical model, (b) A7 vs. Jo. 5 under two 
different experimental conditions. The filled circles are data 
with elevated shot-noise temperature T* , whereas the open 
circles are with variable bath temperature T. The solid line 
is the result of thermal activation model with constant energy 
gap, and the dashed line takes into account the suppression 
of the BCS gap p3| in the measurement with increased T. 



to temperatures T* ~ 5 K, about four times above the 
critical temperature Tq of aluminium, (ii) At all currents 
employed we are well above the threshold To of thermal 
activation: for this particular sample To — 0.3 K. 

Figure |21 (a) shows the width of the switching thresh- 
old, AI, of sample A. Again the experimental (circles) 
and theoretical (line) results are in good agreement be- 
tween each other using the same fit parameters as in 
Fig. [5] In Fig. [3] (b) data on switching of Sample A has 
been plotted in two different experimental conditions, in 
each case as AI vs. Jo. 5- The open circles are data where 
the bath (lattice) temperature has been varied by heat- 
ing the sample stage to several temperatures T in the 
range 0.03 - IK. The full circles are data taken at the 
base temperature of the cryostat (T ~ 30 mK), but by 
varying the equivalent temperature T* injecting different 
levels of Jn . The line closely following the latter data is 
that originating from the thermal activation model. At 
the first sight it is surprising that the truly thermal data 
(open circles) fall far below this line and the other set of 
data. This is, however, accounted for by the fact that 
when increasing the bath temperature close to Tq of alu- 
minium (~ 1.2 K), the BCS energy gap is progressively 
diminishing, thus leading to decrease of the critical cur- 
rent. The line following closely this data set is obtained 
by using the very same thermal activation model, but 
by taking into account suppression of the BCS gap due 
to temperature T 14]. This figure thus demonstrates 
that our thermal activation model is in good agreement 
with the shot noise data, and that it is indeed possible, 
by promoting shot noise, to study Josephson dynamics 
at super- T3 temperatures without suppressing supercon- 
ductivity. 



Finally there are a few more topics to address. We 
found out that the adiabatic models of rocking slowly the 
current bias of the J J do not account for our observations 
on shot noise: the adiabatic models, involving no exci- 
tations, predict exponential dependence in variance-of- 
the-current of the tunnelling rate from the ground state 
through the barrier 0, [n|. Our experimental results 
demonstrate, however, much weaker dependence, which 
we derived in this Letter. The (nearly) resonant model 
works in this case because the circuit responding to the 
shot noise resonates at frequencies, which are of the same 
order as the plasma frequency of the junction. Further- 
more, the low frequency noise does not run through the 
detector junction but rather leaks through the bias line. 
Note that at DC we suppress the detector current due 
to In totally. The second issue is the assessment of an 
MQT detector as an absolute noise detector. Based on 
the invariance of the obtained results, especially in terms 
of results on different scatterers on the same sample, see 
Fig. 2 (a), we suggest that this kind of a detector could 
be made into an absolute on-chip detector of Fano- factors 
0, and noise in general, by careful tailoring of the circuit 
surrounding the JJ and by measuring Q independently. 
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